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1. Introduction 

The study of the boundary behavior of nonnegative solutions of elliptic or parabolic equations 
occupies a central position in the theory of second order pde's. One fundamental result states 
that if it and v are two nonnegative harmonic functions in a Lipschitz domain, and it and v 
vanish continuously on a portion of the boundary, then they must converge to zero at the same 
rate. A quantitative, scale-invariant version of this statement is known as comparison theorem, 
and it was first proved by Dahlberg [D] and Wu [WuJ, and subsequently generalized by Jerison 
and Kenig to non-tangentially accessible (NTA) domains in [JKj. At the same time, in [CFMS 
Caffarelli, Fabes, Mortola and Salsa extended this principle to solutions of divergence form 
elliptic equations in Lipschitz domains, whereas in [BJ Bauman proved the comparison principle 
for solutions of nondivergence form elliptic equations in Lipschitz domains. The proof of the 
comparison theorem hinges crucially on the so-called boundary Harnack principle, or Carleson 
estimate, which states that for a nonnegative harmonic function which vanishes continuously on 
a portion of the boundary, then its values in an interior region and up to that boundary portion 
are controlled by the value of the function itself at a suitably scaled interior point. 

The first version of the boundary Harnack principle for parabolic equations appeared in Kem- 
per's paper |Ke| on the non-tangential convergence of solutions of the heat equation in Lip(l, i) 
domains in R n+ . Such result was subsequently extended by Salsa in [S] to divergence form par- 
abolic equations in Lipschitz cylinders, and by the second named author in [G] to non-divergence 
form parabolic equations in Lipschitz cylinders. In the same paper a parabolic version of the 
comparison theorem was first established for C 2 cylinders. In [G] it was also proved for the 
first-time that nonnegative solutions vanishing on the lateral boundary of non-divergence form 
parabolic equations with time-independent coefficients satisfy a basic elliptic type Harnack in- 
equality, the so-called backward Harnack inequality. In [FGST] the parabolic comparison theorem 
was obtained for the heat equation in Lip(l, ^) domains in M n+1 , whereas in |FGS2j this result, 
the ensuing doubling property of the caloric measure, and the backward Harnack inequality 
for equations with time-independent coefficients were established for divergence form parabolic 
equations in Lipschitz cylinders. In the paper [FSJ Fabes and Safonov first succeeded in proving 
the backward Harnack inequality for divergence form equations with time-dependent coefficients. 
Subsequently, this result was generalized to nondivergence form parabolic equations by Fabes, 
Safonov and Yuan in [FSYJ. A unified approach to the backward Harnack inequality and the 
comparison theorem, which covers at the same time variational and nonvariational parabolic 
equations, was developed by Safonov and Yuan in the paper |SY| . 

In their recent paper [ASS] Armstrong, Sirakov and Smart have established the comparison 
theorem for solutions of fully nonlinear elliptic equations, in fact more generally for solutions to 
differential inequalities for the Pucci extremal operators in C 2 domains. The authors adapt the 
proof of Theorem 1.3 in [BCNj . see Proposition 2.1 in [ASS]. The comparison theorem has been 
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one of the crucial tools in [ASS] in studying the behavior of viscosity solutions of fully nonlinear 
elliptic equations near singular boundary points. The same result for smooth domains has also 
been employed in the recent interesting paper |AS] to establish the strong unique continuation 
property for solutions to fully nonlinear elliptic equations under some structural assumptions. 

These recent results in fully nonlinear elliptic equations, and the desire of investigating their 
parabolic counterpart following the historical development for linear equations outlined above, 
have provided a first natural motivation for us in studying the boundary behavior of nonnega- 
tive solutions of fully nonlinear parabolic equations. In this paper we consider fully nonlinear 
parabolic equations in M n+1 of the type 

(1.1) F(D 2 u,Du,x,t) = u t , 

where F satisfies the structural conditions in (|2.2p below, along with F(0, 0, x, t) = 0. Equations 
such as (jl.ip are the fully nonlinear counterpart of the non-divergence form linear equation 

n n 

(1.2) aij(x, t)DijU + bim = u t . 

i,j=l i=l 

Viscosity solutions to (jl.ip are a subclass of the (viscosity) solutions to the differential inequal- 
ities in (|2.5p below. Such inequalities involve the parabolic extremal Pucci operators and the 
class of their viscosity solutions will be denoted by the notation S(X, A, a), see Definition 12.21 
and more in general Section [2] below for a detailed discussion of this aspect. 

The primary purpose of the present paper is studying the boundary behavior of nonnegative 
functions in the class S(X,A,a). Since, as we have just observed, viscosity solutions to (jl.ip 
are in this class, it follows that our results provide corresponding statements for nonnegative 
viscosity solutions to (jl.ip as a special case. In this respect, we would like to emphasize that the 
results in the present paper encompass a much larger class of equations, than just fully nonlinear 
uniformly parabolic equations such as (jl.ip . For instance, in Section [5] we show an application 
of our results to the following nonlinear singular diffusion equation 

(1.3) \Du\ 2 - p dw(\Du\ p - 2 Du) = u t , 1< p < oo, 

which has recently been studied in [MPRj . [Do] and |BGj . This is an evolution associated 
with the p-Laplacian that generalizes the motion by mean curvature, which corresponds to the 
case p = 1. Viscosity solutions of the equation (|1.3|) satisfy (in the viscosity sense) the double 
differential inequality (|2.5p below, and therefore belong to the class 5(A, A, a), with 

A = min{l,j> — 1}, A = max{l,p — 1}, a = 0. 

However, the equation (jl.3j) does not lend itself to be cast in the form (|1.1|) above. Since 
the results in this paper pertain functions in the general class 5(A,A, a), it follows that they 
encompass the degenerate evolution (jl.3|) as well. 

Apart from their relevance in the historical development of the subject, and their applications 
to degenerate diffusion equations such as (jl.3p . we believe that our results constitute a first step 
in: 

(i) studying the parabolic counterpart of the results obtained in [ASSJ ; 

(ii) establishing space- like unique continuation for solutions to fully nonlinear parabolic equa- 
tions. We mention in this connection the recent paper [Ba] in which the first named 
author has obtained a partial generalization to parabolic fully nonlinear equations of the 
unique continuation result in |AS| : 

(iii) providing some of the basic tools necessary to the study of free boundary problems for 
fully nonlinear parabolic equations. 

This paper is organized as follows. In Section [2] we introduce the relevant notions and gather 
some known results for fully nonlinear parabolic equations. In Section [3l by employing the 
barrier in [M] and using ideas from [G], we prove the Carleson estimate, or boundary Harnack 
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principle, for nonnegative functions in the class <S(A,A, a) in Lipschitz cylinders, see Theorem 
13.31 below. As an application, we establish an elliptic type Harnack inequality for nonnegative 
solutions vanishing on the lateral boundary, see Theorem 13.41 below. In Section H] we establish 
our main result: the local comparison theorem for nonnegative functions in the class 5(A, A, a) 
in C ' cylinders, see Theorem 14.11 below. We show that the barriers in [G] can be adapted to 
functions in the class S(X, A, a). In this context, we would like to mention that although the 
equation (|1.1|) above is not invariant under translations in space and time, however the class 
S(X, A, a) does possess such invariance. It is this simple yet remarkable fact that allows us 
to establish the backward Harnack inequality for nonnegative functions in 5(A, A, a) vanishing 
on the lateral boundary, see Theorem 14.31 below. In closing, we mention that extending the 
comparison Theorem 14.11 to cylinders with a Lipschitz, or even non-tangentially accessible base, 
remains a challenging open problem. To the best of our knowledge, the comparison theorem 
is not known even in the elliptic case. Such results would have important applications to free 
boundary problems for fully nonlinear equations. We intend to come back to these questions in 
a future study. 

Acknowledgment: We would like to thank Adrzej Swiech for bringing reference |CKSj to our 
attention and for a helpful overview of the parabolic viscosity theory. 

2. Preliminary results 
We consider fully nonlinear parabolic equations in W n+1 of the type 

(2.1) F(D 2 u,Du,x,t) = u t , 

with an F which satisfies F(0, 0, x, t) = 0. Hereafter, the notation S n will indicate the space 
of (real) n x n symmetric matrices M. Given a bounded open set f2 C W 1 , and T > 0, by 
saying that (|2.ip is uniformly parabolic in the cylinder Ct = O X (0, T) we mean that there exist 
< A < A, and a > 0, such that for every N G S n , N > 0, every p,q G M n , and (x, t) G Ct, one 
has 

(2.2) raA||iV|| -a\p-q\ < F(M + N,p, x, t) - F(M,q,x,t) < A\\N\\ + a\p - q\, 

where we have let ||M|| = sup|Mx|. 

|z|=i 

Whenever convenient, we will use the notation z = (x,t),Zo = (xo,to),C = (£> r ) to indicate 
points in R n+1 . We now introduce the relevant notion of subsolutions and supersolutions, see 
p. 34 in [WJ and p. 2004 in [CKS]. Let f2 be an arbitrary open set in By a parabolic 

neighborhood of a point zq = (xo,to), we mean the interesection of an Euclidean neighborhood 
U of zq with W 1 x (—oo, to]- I n the following definition, a local extremum is to be understood 
with respect to parabolic neighborhoods. 

Definition 2.1. Let D C M n+1 be a bounded open set. A function u £ C(D) is said to be a 
viscosity subsolution to (|2.ip in D if, for given <p G C 2 (D), at a local maximum zq of u — cp, we 
have 

tpt(zo) < F(D 2 (p,Dp,u,z ). 

A function u G C(D) is said to be a viscosity supersolution to (|2.ip if, for given (p G C 2 (D), at 
a local minimum zq of u — ip, we have 

(2.3) <pt(zo)>F(rP<p,D<p,u,Zo) 

If u is both a viscosity subsolution and supersolution, then we say that it is a viscosity solution 
to d^U in D. 
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Given < A < A, we will denote by 

[[X,A]] = {A€S n \M n <A<M n }. 

We let V^ A and V^ A denote the maximal and minimal Pucci extremal operators corresponding 
to A, A, i.e., for every M 6 S n we have 

V+ A (M) = V+ A (M, A, A) = A e t + A £ e u 

ei>0 e;<0 

P A - A (M) = P" A (M, A, A) = A e * + A E 

e;>0 e;<0 

where ei = e,(M), i = 1, ...,n, indicate the eigenvalues of M. As is well known, V^ A and V7 A 
are uniformly elliptic fully nonlinear operators, and moreover 

(2.4) V+ A (M) = sup [trace(AM)] , V^ A (M) = inf [trace(AM)] . 

Ae[[A,A]] ^e[[A,A]] 

What is important for us, see Proposition 3.10 in |Wj . is that a viscosity solution u to ()2.ip 
satisfies the following differential inequalities in the viscosity sense 

(2.5) V^ A (D 2 u) + a\Du\ - u t > > V^ A {D 2 u) - a\Du\ - u t , 



where a > is the number appearing in (|2.2p . Henceforth, we write 

(2.6) £ + u = V+ A (D 2 u) + a\Du\, C~u = V+ A (D 2 u) - a\Du\. 

Definition 2.2. Given a bounded open set D C the symbol S(A, A, a) will indicate the 

class of functions u G C{D) which are viscosity solutions of (|2.5p in the sense that u is at the 
same time a viscosity subsolution of the fully nonlinear equation 

(2.7) C + u - u t = 0, 
and a viscosity supersolution of 

(2.8) £~u-ut = 0. 

Remark 2.3. We note that, unlike that of solutions to (|2,ip above, the class S(X, A, a) in 
is invariant under orthogonal transformations in the space variable and also, importantly, it is 
invariant under the time translations t — > u(-,t + h), fixed. This simple, yet remarkable, 

fact will play an important role in Section [7] below. 

For further details on the viscosity theory for functions in the class 5(A, A, a) we refer the 
reader to section 3 in |Wj or section 1 in [CKSj . We have observed above that if u is a viscosity 
solution of (|2.ip . then u £ <S(A, A, a) (the opposite implication is not true, see the discussion in 
the introduction, and Section [5.21 below). Since all the results in the present paper pertain to 
functions u £ <S(A, A, a), we infer that our results encompass viscosity solutions of the equation 
flUD as well. 

For equations of the form (|2.ip , we have the following comparison principle (see for instance the 
paper [GGISj . which covers more general degenerate equations and even unbounded domains). 

Remark 2.4. We note that, since the fully nonlinear equations (|2.7p and (|2.8p are special 
instances of (|2.ip . Theorem \2.5\ and Corollaru \2.6\ below are also valid for the operators (|2,7p . 



Theorem 2.5. Let u be a viscosity subsolution and v be a viscosity supersolution in Ct of the 
fully nonlinear equation (|2.ip . If u < v on 8 p Ct, then u < v in Cx . 

The above comparison theorem has the following corollary, which will be subsequently used 
in this paper. 
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Corollary 2.6. Let D C W l+1 be a bounded open set, let u be a viscosity subsolution and v be 
a viscosity supersolution in D of (|2.ip . If u < v on dD, then u <v in D. 

Proof. When D is a parabolic cylinder Ct, then the conclusion follows from Theorem 12.51 above , 
It easily follows that Corollary 12,61 continues to hold for a finite union of such cylinders. In 
order to see this, consider K = |jf=i Qi x (*ii>^2i) and assume that n<uon dK. Take points 
s\ < S2 < S3 < ... < sjy in the increasing order such that s^ = Uj for some i = 1, 2,j = 1, .., M. 
Consider the set K n [s\ < t < S2] which is a disjoint union of sets of the form 0^ x (si,S2) 
where each 0,^ is a connected domain (union of subcollection of cubes Q^s). Moreover, u < v 
on <9 p (f2fc x (s\,S2)) for each k. By Theorem (|2.5[) . we thus have u < v in Qk x (si,^)- We 
conclude that u < v in K n [t < S2]- Now, we proceed inductively to obtain the same conclusion 
in K n (sj <t< Sj+i], and the claim follows. For a general domain L>, given e > 0, consider the 
set 

= G n I u(x, t) > v(x, t) + e}. 

Suppose K 7^ 0. In view of the assumption u < u on 9D, the set X is compactly contained in 
D. Therefore, there is an open set D C D such that K <Z D and L> is a union of finitely many 
cylinders. On dD, we have u(x, t) < v(x, t) + e (note that t> + e is also a supersolution of (|2.7p . 
or (|2.8|) ). From what was said above we conclude that u{x,t) < v{x,t) + e in D, a contradiction 
to the fact that K C D. Therefore, K = 0, and we thus have u < v + e for all e > 0. The 
desired conclusion follows by letting e — > 0. 

□ 

Remark 2.7. Hereafter in this paper when we say that a constant is universal we mean that it 
depends only on the dimension n, and the parameters A, A, a in (|2.5p . 



Now fix i? G (0, 7r), and r > 0. We will denote by Y$ r the truncated cone of radius r, vertex 
at 0, aperture # and axis along the x n coordinate axis, defined by 

T^ r = B r (0) n fix G M n I arccos fej < i?J U {0}^) . 

We denote by r^ j00 = |x G M n \ arccos (^n < j U {0} the infinite cone corresponding to the 

case r = 00. By translating and rotating r$ jr we obtain a truncated cone r$ jr (xo) with vertex 
at any other point xq G W 1 . We need the following barriers on T^ jr for the operator C + . The 
existence of such barriers follows directly from a result of Miller in [M] for non-divergence form 
elliptic operators with bounded measurable coefficients. 

Proposition 2.8. Given G (0,7r), i/iere exist a universal Ro > 0, depending also on and 
a uniform barrier ip for r# ^ and the operator C + . This means that <p G C 2 in the interior of 
F^,_r and that we have for some universal C,a,fj,i,K > 0, depending also on &: 

1) < tp(x) < C\x\ a ; 

2) ^(0) = 0; 

3) ip(x) > m, \x\ = R ; 

4) £+<f(x) < -K\x\ a ~ 2 . 

Proof. The barrier in Theorem 2 in [M] is of the form (p(x) = \x\ a h ^arccos It is proved 

in [M] that, for all matrices [a^] G [[A, A]], one has in r^^ \ {0} 

OijDijipfa) < -A\x\ a ~ 2 , \Dip(x)\ < -Blx)"" 1 , 

where A, B > are universal constants depending also on Since by fj2.4j) we have 

^A^ 2 ^ 9 )^) = SU P aijDijip(x), 

Xl n <[aij]<Al n 
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we thus obtain from (12.61) 



C+ip = V+ A (D 2 ip) + a\D<p\ < -A\x\ a ~ 2 + aB^- 1 = -^|x| Q " 2 [l - A _1 oB|x|]. 

Since |x| < Rq for x £ ^-&,r , it suffices to choose for instance Rq < ^° obtain the desired 
conclusion with K = A/2. 

□ 

Now, if ry r indicates a truncated cone of radius r, aperture #, and vertex at 0, we denote by 

cv = iv x (-^V 2 ) 

the parabolic wedge with vertex at (0,0) £ radius r, aperture i? and height 2r 2 . Us- 

ing Proposition 12.81 and adapting Theorem 1.12 in [G] for linear parabolic operators in non- 
divergence form, we can construct a barrier for C$ r and the fully nonlinear parabolic operator 
£ + -d t . 

Theorem 2.9. There is a universal Rq < 1, depending also on 9, such that for < r < Rq, 
there exists a uniform barrier for C + — dt on C# r . This means that we can find a function 
tp £ C 2 in the interior of C$^ r for which we have for some universal C,a,fi > 0, depending also 
on i9: 

1) o<^x,t)<c(^f^y ; 

2) 0(0,0) 0: 

3) tp(x,t) > \x for each (x,t) £ dC$ r such that either \x\ = r, or t = — r 2 ; 

4) £+(V) - i> t < 0. 

Proof. We begin by observing that the theorem holds true when r = Rq where Rq is as in 
Proposition (|2.8|) and then analyze the general case. Without loss of generality, we can assume 
R < 1. We define 

K 

iP{x,t) = tp{x) + —\t\ K , 

where k£N with k > 2. From Proposition 12.81 above, by arguing as in the proof of Theorem 
1.12 in [G], it is easily seen that we can reach the desired conclusion for r = Rq. For the general 

R R? 

case, we define ipi(x,t) = tp(-f-x, for (x,t) £ CV, r - Let us notice that (x,t) £ C$ }r if and 

R i? 2 

only if (-^x, ^t) £ CV,_r . Then, we have 



V+ K {D 2 ^) + — a\Dij}\ - i> t 



■Rq R%} 
„2 L J- 



— X, °' 



i?0 

R R 2 

Now, when < r < Rq we have < 1, and since (yi, ^-t) £ (7$^, by the properties 4) of ^ 
in Ci9 5 ij we conclude 

[£+(^) - (^) t ] (s,t) < 5f [V+ A (D 2 1P) + a\Di/>\ - ^] {^x, ^t) < 0. 

The properties l)-3) are also easily verified for the function ip\ in the wedge C$ jr . Therefore, 
the conclusion follows. 

□ 

We close this section by quoting the important Krylov-Safonov type Harnack inequality for 
solutions of the differential inequalities ([23]). For this result we refer the reader to Theorem 
4.18 in |Wj . when no lower order terms are present, and to Section 4.6 in the same paper for a 
discussion on how to handle lower order terms. One should also see Lemma 5.2 in |CKS] for the 
explicit statement. For rj > 1 and r £ [0, 1], define 

Q(rj,r) = {(x,t) £ R n+1 | max \xi\ < r, < t < r?r 2 }. 

Ki<n 
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Theorem 2.10. Let u > be a function in the class S(X, A, a) in Q(rj, r) and < a < 1. Then, 
there exists a universal constant C = C(r],a) > such that 

(2.9) max u(x,r 2 ) < C min u(x,i]r 2 ). 

\x\<crr \x\<crr 

3. Boundary Harnack Principle for Lipschitz cylinders 

In this section, we work on a cylinder Ct = Ox (0,T), where O C W 1 is a Lipschitz domain. 
We recall that a bounded domain VL is a Lipschitz domain if for each Q G dQ there exist a ball 
B TQ (Q) centered at Q and a coordinate system for R n such that in these coordinates 

B ro (Q)nn = B ro n{(x',x n ) ■. x' eR n ~\x n xt>(x')}, 

B ro (Q)ndn = B ro n{(x',x n ) : x' G K n ~\x n = <j)(x')}, 

where <f> : M ra_1 — )• 1R is a Lipschitz continuous function with 1 1-0011^(^-1) < m. By compact- 
ness we can, and will, assume that the radius ro and the constant m are independent of Q G dtl. 
Such constants determine what will be called the Lipschitz character of Q. We further indicate 
with S T = dflx (0, T) the lateral boundary of C T , while d p C T = S T U (ft x {0}). For (Q, s) £ S T 
and sufficiently small r < j^, we define 

^ r (Q,s) = {(z,i) G C T I \x - Q\ < r, \t - s\ < r 2 }, 

A r (Q,s) = y r (Q,s)r\d p C T . 
If Q G dtl is given by (x' Q , ^(x' )) in the above local coordinates, for r > small enough we set 

A r (Q, s) = (x' , <t>(x' ) + r,s + 2r 2 ), A r {Q, s) = (x' , <f>(x' ) + r,s- 2r 2 ). 

A Lipschitz domain 0, has a uniform exterior (and interior) cone T$ t00 (Q) at every point 
Q G d£l. The uniform aperture •& G (0,7r) of such cones depends on the Lipschitz constant m of 
the local defining function </> of dfl. 

By adapting the proof of Lemma 2.1 in [Gcj . we next show that a solution to the fully nonlinear 
equation (|2.5|) which continuously vanishes in a surface neighborhood of a point on St is actually 
Holder continuous up to the boundary. 

Lemma 3.1. Let (Q,s) G St- There is a universal R$ > 0, depending also on the Lipschitz 
character of Q, such that for < r < Rq, and sufficiently small depending also on s,T — s, ro, 
and any u > in the class S(X,A,a) in *$> r (Q,s) and vanishing continuously on A r (Q,s), we 
have for universal C,a > 0, depending also on the Lipschitz character ofQ,, 

(3.1) u ( I)f )<c|tMl^| M ^ (ar,t)€* r (Q,«), 

where M r (u) = sup u. 

*r(Q,s) 

Proof, as in Proposition 12.81 

Without loss of generality, we may assume that M r {u) < oo. Let T$ yOQ (Q) be a fixed infinite 
cone exterior to O and having vertex at Q. If T r / 2 (Q) = B r /2(Q) \ r^ i00 (Q), let C r / 2 (Q,s) = 
F r / 2 (Q) x (s — r 2 , s + r 2 ) denote the parabolic wedge at (Q, s), and let be a uniform barrier 
on C r /2(Q, s) constructed as in Theorem 12.91 above. In the cylinder Ct H C r / 2 (Q, s), we consider 

. . M r (n) ,. . 
w[x,t) = ip\x,t), 
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where \i is as in Theorem 12.91 By 4) in that theorem we have C + (w) — wt < in CtH C r / 2 {Q, s), 
whereas 3) implies w > u on 8 p (Ct H C r / 2 (Q, s)). Since u is a subsolution and w a supersolution 
of C + — dt , by the comparison principle Theorem 12.51 see Remark 12.41 we conclude 

(3.2) u(x,t) < Mr ^ ip{x,t), (x,t) eC T nC r/2 (Q,s). 

From 1) in Theorem 12.91 we obtain the desired bound (pTXJ on u in Ct H C r /2(Q, s). If instead 
(x,t) G ^ r (Q,s)\C r/2 (Q,s), then \x-Q\ > r/2. Hence, in that set , we have ^'-Ql+I^l 1 / 2 ^ > 
2~ a . Therefore for (x,t) G Vl/ r (Q, s) \ C r / 2 (Q,s) we trivially have 

u(x,t) < M r (n) < 2 a ^ |x ~ 9l+ ^ |t ~ s|1/2 ^ M r (u). 
This completes the proof. 

□ 

The next statement provides a quantitative information on the growth near the boundary of 
a nonnegative solution. We will use this result to start the worsening process in the proof of 
Carleson's estimate, Theorem 13.31 below. Its proof is based on the non-tangentially accessible 
character of a Lipschitz domain, and on the Harnack inequality in Theorem 12.101 above. Once 
these tools are available, one can carry the argument similarly to the proof of Lemma 2.2 in [C-] , 
We thus omit it, and refer the reader to that source for the relevant details. 

Lemma 3.2. Let (Qo,so) G St, ol G (0,1). There exists a universal constant C = C(m,a) > 
such that, if u > is a function in the class <S(A,A, a) in ^f2r(Qo,so), then the following 
condition holds: if (y,s),(x,t) G ^ 3r / 2 (Qo, sq) are such that d\st(x,dD) = ar, t — s > r 2 /2, 
dist(y, 3D) > where h G NU {0}) is such that 2~ h < a, then 

(3.3) u(y,s) < C h u(x,t). 

With Lemmas 13.11 and 13.21 in our hands, we can argue as in the proof of Theorem 2.3 in [G] 
and establish the following Carleson estimate, or boundary Harnack principle. 

Theorem 3.3 (Carleson estimate for Lipschitz cylinders). Let (Qoi s o) G St- Then, there exists 
a universal constant C = C{m) > such that, if u > is a function in the class S(X, A, a) in 
^2r(Qo)So) which vanishes continuously on A2 r (Qo) so), then 

(3.4) u(x, t) < Cu(A r (Q , s )), {x, t) G ^ r/8 (Qo, «o), 
where r is small enough depending on sq,T — So,7"o- 



Proof. We argue as in [U] and we first conclude (|3.4p for r < Rq where i?o is as m Lemma (|3.ip . 
Consequently, by an application of interior Harnack inequality in Theorem 12.101 we conclude 
that validity of (|3.4p for r depending only on ro, T — so, sq for a possibly larger universal constant 
C > 0, depending also on the Lipschitz character of O. 

□ 

As an application of the Carleson estimate, we have the following elliptic type Harnack in- 
equality when a solution to (|2.5p vanishes continuously on the lateral boundary St of the cylinder 
Ct- 

Theorem 3.4. Let u > be a function in the class S(X, A, a) in Ct- For 5 > set = {x G 
D | dist(z, dfl) > 5}, and C$ : t = ^jX (^ 2 jT). There exists a universal constant C > 0, which 
also depends on m, 5, diam £l,T, such that, if u vanishes continuously on St, then 

(3.5) max u < Cminn 
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Proof. We note that the assumption guarantees that u G C(Qs,t), and that u is a subsolution 
to C + — dt and a supersolution to £~ — dt in O^r respectively. Let 

ki = min it, and k<i = max it. 

d p fls,T d p Qs,T 

Clearly, k\ is a solution (hence a subsolution) to C~ — dt and ki is a solution (hence a superso- 
lution) to £ + — dt- Since k\ < u and ^ > u on d p Qs,T, by the comparison principle, Theorem 
(|2.5|) . we obtain 

ki = min it, = max u. 

Moreover, by the compactness of d p Qs,T, there exist (Xq,Tq), (X\,Ti) £ d p £ls,T such that 
u(Xq,Tq) = min u = min it, and u(X\,T\) = max u = maxu. 

The rest of the proof now follows from Theorem 13.31 and the interior Harnack inequality in 
Theorem 12. 10|. as in the proof of Theorem 2.6 in [G]. 

□ 

4. Comparison theorem for non-negative solutions in a C 1,1 cylinder 

In this section we establish a comparison theorem for functions in the class <S(A,A, a) in 
Ct = fix (0,T), where fi C M ra is a C ' domain. Our main result, Theorem Boundary 
comparison theorem, should be compared with Proposition 2.1 in [ASS] in which the authors 
establish its elliptic counterpart for C 2 domains. 

We recall that a C ' domain is characterized by the existence of a uniform exterior and 
interior tangent ball. The largest radius of such balls determines what we call the C 1 ' character 
of the domain fi. Our proof of the comparison theorem will follow by adapting the barriers in 
[G] to the fully nonlinear case, and arguing similarly relying on the natural sublinearity, in place 
of linearity. The following is our main result. 

Theorem 4.1 (Local comparison theorem for C 1 ' cylinders). Let (Qq,sq) £ St- Then, there 
exists a universal ro > 0, depending also on the C 1 ' 1 nature of SI and on sq,T — so, and a 
universal constant C > 0, depending also on the C ' character oftt, such that if for < r < r^, 
u, v > belong to 5(A,A, a) in ^2r(Qo,so), and u = v = continuously on A2 r (Qo,so), then 
for all (x,t) 6 ^ r (Qo,so), one has 

, 4] x 1 u(A r (Q ,s )) < u(x,t) < c u(A r (Q ,s Q )) 

C v(A r {Q ,s )) ~ v(x,t) ~ v(A r (Q ,s ))' 

Proof. For Q G dQ, let vq be the interior unit normal to d£l at Q and for r > set 

(4.2) ^Q) = Q + L UQ £ 2 (Q) = Q- Yq u Q- 
Let us define for 5 > 

Ps(x,t) = {(y,t) | \x-y\ 2 + \t-r\ < 5 2 }. 

We now choose and fix ro > 0, depending on the C > character of Q,, so that for < r < ro, 
and each (Q, s) E A r (Qo, so) we have 

(4.3) < v Q , u Qo >> -, for Q G dn, \Q - Q \ < r, 
and 

(4.4) p r _ (^(Q), s )ns T = Ps. (6(g), s)ns T = {(Q, s)}. 

From now on, (Q, s) is assumed to be a arbitrary point of A r / 8 (Qo) so)- We consider the set 

5 r (g,fl) = P i (Q, a )ni' i (ei(Q),«). 
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Adapting the idea in [G] . we wish to find two suitable barriers h and / such that 

C~{h) -h t >0 in S r (Q,s), 

and 

£ + (f)-ft<0 in MQo,s )nP r/s (UQ), s ). 

With these inequalities in place, by an application of the comparison theorem (Corollary 
above), we then conclude the following estimates 

u 

(4.5) — — - >h in S r (Q,s), 

u{Ar(Q0,S0)) 

and 

(4-6) < Cf in * r n P r/8 (Q, s). 

u{A r (Qo,s )) 



We emphasize that (|4.6p is the crucial step. To accomplish it we will use the Carleson estimate 
(|3.4p in Theorem 13.31 above. 

We assume henceforth that r > is chosen sufficiently small (depending on tq, sq and T — so) 
so that all the forthcoming arguments make sense. Observe that if 

ti = inf{t | (x, t) E S r (Q, s)}, t 2 = sup{t | (x, t) E S r (Q, s)}, 

then there exist a±, a 2 > such that 

(4.7) s - 2r 2 + air 2 < ti < t 2 < s + 2r 2 - a 2 r 2 . 
We next consider the sets 

$ 1 r (Q,s) = dS r (Q,s)ndPz(( 1 (Q),s), <Z> 2 r (Q,s) = dS r (Q,s)ndPr_(Q,s). 
There exists 03 > such that 

(4.8) dist($ 2 , St) > a 3 r. 

Therefore, from (|4.7p . (|4.8|) and the Harnack inequality in Theorem 12.101 ab ove . we obtain 

(4.9) u(x, t) > Cu{A r (Q , s )), (x, t) E $ 2 (Q, s). 

Now, let R(x,t) 2 d = |x-^i(Q)| 2 + \t-s\. Observe that on &P r / 4 (£i(Q), s) we have R 2 = r 2 /16, 
and that 

(4.10 mm = = — 

S r (Q,s) r 2 64 

Now we choose a > large enough depending on A, A, n. Let 

/ / aR 2 \ / a 
Gn = max exp jr- — exp 

We stress that, in view of (|4.8p . Co is independent of r > 0. Pick 7 > such that 

(4.11) ^<^' 

2Go 

where C > is the constant in (|4.9p above. We now define 

aR 2 \ ( a 

— - eXP (-16 



/i(x, t) = 7 ( exp ( — ) - exp 



A verification shows that h = on &l(Q, s). Moreover, from the choice of 7 and ()4.9|) . we have 
that 

(4.12) MM)<7Q)<f <<?< , v (x,t)E$ 2 (Q, S ). 
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Now, for (x,t) £ S r (Q,s) we have 

(4.13) Dh = -^(x-£ 1 (Q))exp(-^- 

and 

aBf 



(4.14) Dijh = ( Xi - ei(Q)i) ( Xj - UQ)j) ~ ^T<%) exp (- 



In view of (|2.4p above we thus find 
(4.15) 



2ct"f ( 2a n \ ( aR^ \ 

Therefore, from (|4.10p . and by choosing a large enough (depending also on the C 1 ' character 
of S7), we have in S r (Q, s) 



(4.16) V- XA (h) > 5 (^\x - e l( Q)| 2 - nA) exp (-^ 



l V r 



2cry / 
> — 7T- exp 



Hence, from (|2.6|) and (|4.16p we find in S r (Q,s) 

(4.17) £-(/0>^(l-a|z-£i(Q)|)exp / ' ni? " 

Now, if (x,i) G S r (Q,s) we have |rr — £i(Q)| < r, and therefore if we choose < r < i we 
obtain 

MS) ^C)£|?exp(-^ 
On the other hand, we have 

cry . . / ai? 
Aii = — — sgn(t - s) exp 

Therefore, 



(4.19) - /i t > ^Jexp (-^r) (3/2 - 1) > 0, 

which shows that h is a subsolution to £ _ — dt in the set (t ^ s), i.e., almost everywhere in 
S r (Q, s). Hence, h is a strong subsolution in S r (Q, s). By Proposition 2.11 in [CKSJ, we conclude 
that h is also a viscosity subsolution of £~ — dt in S r (Q,s). Since h(x,t) < gjfe^ji on <9SV, 
and by the assumption that it £ <S(A,A, a), u is a supersolution of the same operator, by the 
comparison principle Corollary 12.61 applied to C~ — dt, we obtain that 

(4.20) h(x,t)< U( - X n t] in S r (Q,s). 

u(A r (Qo,s )) 

At this point, by the argument on p. 290 in [G] . which uses also (|4.3j) above, we conclude that 
there exists C\ > 0, depending on n, A, A and the C ' character of <90, such that 

(4.21) h(x,s) > c j*"^ 
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when x = Q + Srv Qo ,5 G [0, 1/32], see (3.20) in [GJ. From (PU|) . (IOT|l . it follows that for such 

x 

(4.22) Cl l -?—^u(A r (Q ,s )) < u(x,s). 

r 

The estimate (|4.22p thus establishes the linear growth from below near the lateral boundary St 
for u. 

Our next objective is proving a similar estimate, but from above. For this, we consider the 
set ^> r (Qo, so) n P r /8(&(<9), s), and define 

Ri(x,t) 2 = \x - &(Q)\ 2 + \t - s\, 

where &(Q) is defined in (|4.2p above. If (x,t) S ^ r (QoiSo) H Pr/si&iQ), s), we notice that for 
some constant C5 > we have 

(a o*\ k-6(Q)l 2 > r 

(423) i^M) 2 ^ Cs - 

For k > to be subsequently chosen large enough, depending on n, A, A, C5, we define 

/(M)=1 - 7 Gdb)*) ' 

with 7 = (2gg) fc - Since R\(Q,s) 2 = this choice of 7 ensures that f{Q,s) = 0. In the next 
formulas, to abbreviate the notation we write R\ instead of R\(x,t). A computations gives 



Df = 2jk 



r 2 \ k x-UQ) 



K>2 / r>2 



(2 \ fe 7^ , -] / 2 \ ^ x 

Therefore, (|2.4p above gives 

(4.24) p+ A ( / ) = 2 ^^y sup /_ 2 ii+l) aij ( x _ 6 (Q) )j ( x _ 6 (Q) )j+aii ] 

^1 \Ri/ Ae[[A,A]] I R\ J 

From (|4.24p we easily obtain 

v U f) -%^\) {-^r Ak " 6(g)|2 + nA l • 

Because of (|4.23p . we can now choose k large enough so that, for (x, t) E ^ r (Qo, so)nP r /g(^2(Q), s), 
we have 



R^ \ R 

This estimate and (|2.6p above give 

(4-25) £+(/)< -^(^) (l-a|x- 6(0)|). 

Observe now that, for (x, t) £ \tv(Qo> so) H Pr/si&iQ), s), we have |x — 6(Q)| < i?i < r. Hence, 
if < r < ^, we obtain from (|4.25p 

£+(/)< "^(^) > in* r (Qo,*o)np r/8 (6(g), a ). 



BOUNDARY HARNACK ETC 13 

Since 

ft = ~R% \¥j sgn(t ~ s) ' 

we conclude that C + (f) — ft < when (t ^ s) in ^ r (Qo, sq) n -P r /sfe(Q); s). Thus, / is a strong 
supersolution, and hence a viscosity supersolution of C + {f) — ft, by Proposition 2.11 in |CKS| . 
Moreover, on ^ r (Qo) so) H dP r / 8 ((,2(Q), s), we have 

(4.26) /(z,i) = l-2- 2fc , 

while on A r (QojSo) HPr/g^fQ)^), we have f(x,t) > (as (Q,s) is the closest point where 
it is zero). The crucial step is to now invoke the Carleson estimate, Theorem 13.31 above, which 
gives 

(4.27) < C, (x,t) g MQo,s Q )nP r/s fa(Q),s). 
u(A r (Qo,s ) 

Thus, if C = C(l - 2- 2fe )" 1 , then < C"/ on 9 (* r (Q ,s Q ) n P r/8 (&(Q),s)). Since 

from what we have proved above the function C'f is a viscosity supersolution of C + — dt in 
^r(Qo,so) n Pr/si&iQ), s), whereas, by the hypothesis that u G <S(A, A, a) in Vl/2r(Qo> sq), the 
function ,-r " is a subsolution of the same operator in that set, by C or pilar v 1 2 . 6 1 we conclude 

u{A r {Qo,s ) * 

that 

(4.28) -=r— < C'f in V r (Q , s ) n P/ 8 (6(Q), a). 

u(A-(Qo,so) 

The estimate (|4.28|K combined with the equation (3.33) in [G) . shows that 

b — 01 — 

(4.29) u(x, s) < cV ^u(A r (Q , s )), 

r 

when x = Q + 5ruQ ,5 G [0, 1/32]. Combining (|4.22p and (|4.29p we have proved that 

(4.30) Cl ^^-u(A r (Q ,s )) < u(x,s) < C 2 l -^^u(A r (Qo, so)), 

when x = Q + 5rvQ ,5 G [0, 1/32]. In other words ()4.30p holds all over a slice of ^(Qo^o) 
which is attached to the lateral boundary St, and goes inside Ct for a depth proportional to r. 
The important information (|4,30p tells us that, in a C ' cylinder, any nonnegative function in 
the class <S(A, A, a) in ^2r(Qo,to), which vanishes continuously on A2 r (Qo, sq), must vanish at 
precisely a linear rate. 

We next apply (|4.3(jp to the function v, obtaining 

r 1 1 1 
(4-31) C 2] ~~ — = < < Ci 



x - Q\v(A r (Q ,s )) v(x,s) \x - Q\ vfAriQo, So))' 

when x = Q + 5ruQ ,5 G [0, 1/32]. The estimates (|4.30p . (|4.3ip immediately imply the desired 
conclusion (|4.ip for points (x,s) G ty r (Qo,ta) such that x = Q + 5ruQ , 5 G [0, 1/32]. This allows 
us to now move away from the lateral boundary St by a distance proportional to r. Once this 
is accomplished we can apply the interior Harnack inequality in Theorem I2.1UI above, and thus 
obtain (|4.1|) in the remaining portion of the set ^f r (Qo,to). Hence, we obtain (|4.1|) in ^ r (Qo,to) 
when r additionally satisfies the assumption r < t-. Consequently, by an application of the 
interior Harnack inequality in Theorem 12.101 to both u and v, we conclude that (|4.ip continues 
to be valid for all r small enough depending only on C ' character of f2, sq,T — sq and for a 
possibly larger universal constant C > 0. 

□ 
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As an application of Theorems 13.41 and !4. 11 we obtain the following global comparison theorem 
for nonnegative functions in the class <S(A, A, a) which vanish on the lateral boundary St- This 
result is a consequence of Theorem 13.41 and Theorem 14.11 Since its proof is identical to that of 
Theorem 3.2 in [G], we omit it altogether. 

Theorem 4.2 (Global comparison theorem). Let Ct = O x (0, T) be a C ' cylinder, and 
u,v > be functions in <S(A,A, a) in Ct, such that u,v = continuously on St- Fix a point 
Xq G ft, and for a universal 5 G (0,y/T/2) suitably small (depending on the C 1,1 character of 
Q), set F = ft x (25 2 ,T — 5 2 ). There is a universal constant C > 0, depending also on 5, diam 
ft, and T, such that for all (x, t) £ F one has 

(4 321 1 u(X ,T) < u(x,t) ^ u(X ,T) 

{ ' ' Cv(X ,T) ~ v(x,t) ~ v(X ,T)' 

Before stating the next result we recall that functions in the class 5(A,A, a) are invariant 
under translation in time, see Remark 12.31 above. Consequently, by arguing as in [Gj, Theorem 
3.3, we have the following backward Harnack inequality for functions vanishing on the lateral 
boundary St- 

Theorem 4.3 (Backward Harnack inequality for C 1 ' cylinders). Let Ct, u, Xq £ fl, 5, F and 

C be as in Theorem \4-%\ Then, for r small enough, one has for (x,t) £ F 

(4.33) u(x, t + 4r 2 ) < Cu(x, t). 

Proof. By Remark 12 .31 the function v = u(x, t+4r 2 ) belongs to S(X, A, a). The desired conclusion 
then follows from Theorems ()4.2p and ()3.4p . 

□ 

Remark 4.4. We would like to emphasize that we could have treated slightly more general 
differential inequalities of the form 

(4.34) ^a( u ) ~ a\Du\ - b\u\ - u t < u < V^ A (u) + a\Du\ + b\u\ - u t . 

One can see that the barrier in Theorem \2.9\ can be adapted to this situation as well, and thus 
the Carleson estimate in Theorem \°J.'J\ and the comparison Theorem \4-l\ continue to be valid in 



this situation. However, for the backward Harnack inequality Theorem we crucially make 
use of the boundary maximum principle Theorem \3.4\ which in turn makes use of the fact that 
constants are solutions to the differential inequalities (|2.5p . But this is true in the case when 
6 = 0. Thus, for Theorem \4-3\ we do need the assumption 6 = 0. 



Remark 4.5. Lt remains a challenging open question whether, for functions in the class <S(A, A, a), 
Theorem \4-l\ continues to be true for a cylinder Ct = Q x (0, T), with a Lipschitz, or more in 
general non-tang entially accessible (NTA), base ft. To the best of our knowledge, this result is 
not known even in the case of fully nonlinear elliptic equations. 

5. TWO APPLICATIONS 

In this section we present two applications of the results in Sections [3] and HI We plan to come 
back to further applications in a future study. 

5.1. Non-divergence form linear equations. The spirit of this subsection is mainly demon- 
strative. What we mean by this is that, although the result that we obtain as an application 
of the results in this paper is not new, the way we obtain it is, and puts it in the broader 
perspective of the viscosity fully nonlinear theory. In a cylinder Ct = O X (0, T) consider the 
time-dependent non-divergence, uniformly parabolic linear equation 

n 

(5.1) Pu = Lu — Ut = a ij{ x ^)F>ijU — Ut = 0, 



BOUNDARY HARNACK ETC 



15 



where for (x,t) G Ct one has a,ij(x,t) G [[A, A]] and, in addition, the are assumed to be 
continuous. It is well-known, see for instance Theorem 1.5 in [Gj . that if u is a strong solution 
to d5U, then it G %; ioc (C T ) for all p > 1. Since fl57Lj is a special case of f)2. 1 [) . by Proposition 
2.11 in [CKSJ, u is a viscosity solution to ()5.ip . Consequently, the results in this paper imply 
the backward Harnack inequality, Theorem 14 . 3 1 ab ove . for non-negative solutions to (|5.ip in C ' 
cylinders. We cannot stress enough the non-trivial and remarkable aspect of this observation: we 
are dealing with a linear parabolic equation, ()5.ip . whose coefficients are time-dependent. This 
serious obstruction can be circumvented by looking at (|5.ip from the point of view of the fully 
nonlinear class S(X, A, a) of the Pucci extremal operators, which are instead time- independent! 

Now, having said this, we ought to nonetheless mention that, for the linear equation (|5.ip . 
the backward Harnack inequality in Theorem 14.31 is in fact available, even in Lipschitz cylinders. 
This is so thanks to the work, subsequent to [G], of Fabes, Safonov and Yuan [FSYJ, see also 
the paper [SY] . both of which involve a great deal of hard analysis. In connection with these 
works, we mention that, when u is a viscosity solution of the fully nonlinear equation 

(5.2) F(D 2 u,x,t) -u t = 0, 

where F is concave in D 2 u and C , then, by the Evans-Krylov theorem, u is a classical H 2+a 
solution. This follows, for instance, by the existence of classical solutions with continuous 
boundary values for such F, see Theorem 14.10 in [L], combined with the uniqueness of solutions. 
Therefore, u solves a linear equation such as (|5.1|) . where 

/l Qp 
-^-(sD 2 u(x, t),x, t)ds. 

Hence, from the linear theory, see Theorem 3.7 in [FSYj . in the case when F is concave in M, 
then for viscosity solutions of (|5.2p the Carleson estimate and the backward Harnack inequality 
are true in Lipschitz domains. 



5.2. On a degenerate diffusion equation. In this subsection we show how the results in this 
paper can be applied to obtain new theorems for the nonlinear degenerate diffusion equation 

(5.3) ut = \Du\ 2 - p div(\Du\ p - 2 Du), 1< p < oo. 

The equation f)5.3[) has recently been studied in [BG], |Do| . [MPR| . Such an equation gener- 
alizes the motion of level sets by mean curvature, which is relative to the case p = 1, and the 
heat equation, which corresponds to p = 2. In the interesting paper |MPR] . solutions to the 
equation ()5.3p have been characterized by asymptotic mean value properties. These properties 
are connected with the analysis of tug of war games with noise in which the number of rounds 
in bounded. The value function for these games approximate a solution to the PDE ()5.3p when 
the parameter that controls the size of possible steps goes to zero. Let us introduce the relevant 
notion of sub-, super-, and solution to ()5.3p . For further details on the properties of viscosity 
solutions of (|5.3p . we refer the reader to the paper [BG]. 

Definition 5.1. A function u G C(S1 x [0, T)) n L°°(Q x [0, T)) is called a viscosity subsolution 
of (|5.3p provided that for any ip G C 2 (Q x (0, T)) for which 

(5.4) u — (p has a local maximum at zq G £1 x (0,T], 
then, either 

/ 5 5) JV* < (<% + (p - 2 ) £ 0$ E ) D iW at z o, 

[if Dip(z ) ^ 0, 
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or 

{inf (ipt - (<% + (p - 2)aiaj)Dijip) < at z , 
\a\=l 
ifD<p(z ) = Q. 

Analogous definitions for supersolution, or solution. 

Note that, from Definition 15.11 at a maximum point of u — ip, we have that 

(5.7) ip t < aij(Dip)Dijtp, 
where 

(5.8) min{l,p — 1} I n < [ay] < max{l,p — 1} I n . 

The inequality (|5.7p is seen to be true by considering the matrix [ay] in both cases when Dip 
vanishes and when it does not, sec [BGJ. It is thus easily seen from Definition 15.11 and (|5.8p 
that any viscosity solution to (|5.3p belongs to the class <S(A, A, a), with A = min{l,p — 1} 
and A = max{l,p — 1}. Consequently, and quite notably, all the results in this paper, such 
as the Carleson estimate in Theorem 13.31 the local and global comparison Theorems 14.11 and 
14.21 and the backward Harnack inequality in Theorem 14.31 are valid for nonnegative solutions 
of the degenerate diffusion equation (|5.3p as well. Here, we confine ourselves to state, as an 
example, the following two results. We note that the degenerate parabolic equation (|5.3p has a 
non- variational structure and therefore the viscosity framework of this paper is the natural one 
for studying it. 

Theorem 5.2 (Local comparison theorem). Let Ct = 0, x (0,T) be a C 1,1 cylinder and let 
(QojSo) £ St- For any 1 < p < oo, there exists ro > 0, depending on n,p, the C 1, nature of 
$7 and on sq,T — sq, and a C > 0, depending on n,p and the C 1 ' character of £2, such that 
ifu,v > are two viscosity solutions of (|5,3p in ^?2r(Qo, s o), an d u = v = continuously on 
^■2r{Qo-, so), then for all (x,t) £ ^ r (Qo, sq), one has 

1 u(A r (Qo,s )) < u(x,t) < u(A r (Q ,so)) 
C v(A r (Q , s )) ~ v(x,t) ~ v(A r (Q ,s ))' 

Theorem 5.3 (Backward Harnack inequality). Let Ct, Xq £ Vt, 5, F and be as in Theorem 
\4-^\ above. Then, given a viscosity solution u > of ()5.3[) in Ct, vanishing continuously on St, 
there exist ro > 0, depending on n,p,6, and and the C 1 ' character ofVt, and a constant C > 0, 
depending on n,p and the C 1 ' character of such that for < r < ro, one has for (x,t) G F 

u(x,t + 4r 2 ) < Cu(x,t). 

We close by remarking that, in Definition 15. 11 similarly to what was done in [Do| . [BGJ, [MPR] . 
parabolic Euclidean neighborhoods are considered for local extrema. However, by arguing as in 
Lemma 1.4 in |CKS| (this lemma shows that considering parabolic neighborhoods is equivalent 
to taking Euclidean neighborhoods), it suffices to consider parabolic neighborhoods at Zq. This 
is what is normally done in the literature, see for instance |Wj . 
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